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In  an  earlier  communication  (dealing  with  flow  in  curved  pipes)  a 
numerical  procedure  capable  of  calculating  variable  property  flows  in  arbi¬ 
trary  curvilinear  orthogonal  coordinates  was  reported^.  Convective  and 
diffusive  base  vector  variation  source  terms  appearing  in  the  differential 
momentum  equation  were  tabulated  but  their  derivation  was  not  described.  In 
addition,  and  because  of  the  cases  studied,  only  brief  mention  was  made  with 
regard  to  the  extra  source  terms  arising  in  variable  property  flows.  The 
purpose  of  this  letter  is  to  provide  a  concise  derivation  of  the  differential 

t 

momentum  equation  which  is  the  basis  for  the  numerical  procedure  described  in 
web¬ 
reference  l~and  to  tabulate  general  expressions  for  all  the  sources  which  can 

arise.  From  a  knowledge  of  the  metric  tensor  in  the  curvilinear  system  of 
interest  the  forms  of  these  terms  are  readily  derived  from  the  tables  for  cal¬ 
culation  purposes. 

Throughout  this  letter,  standard  tensor  notation  with  the  summation  con¬ 
vention  is  employed.  A  detailed  exposition  of  the  rules  and  uses  of  tensor 
analysis  in  fluid  mechanics  is  available  in  reference  2.  Briefly,  superscripts 
refer  to  contravariant  components,  subscripts  represent  covariant  components  and 
commas  denote  covariant  differentiation.  Any  quantity  followed  by  an  index 
within  parentheses,  such  as  v(i),  denotes  the  ith  (in  this  case)  physical  com¬ 
ponent  of  that  quantity.  While  the  present  analysis  is  limited  to  steady  state 
forms  of  the  conservation  equations,  it  is  readily  extended  to  time  dependent 
situations. 


t 


2 


The  Momentum  Equation 

The  steady-state  vector  equation  of  motion  for  a  variable  property 
flow  is  given  by: 

In  equation  (1),  ?  is  a  body  force  field  such  as  that  arising  due  to 
gravitational,  buoyancy,  electrostatic  effects,  etc,  ?  is  the  stress  tensor 
field,  and  P  is  the  pressure  field.  The  symbol  denotes  the  divergence 

t 

operator  in  vector  notation,  whereas  $  denotes  the  gradient  operator.  An 
equivalent  tensor  form  in  orthogonal  coordinates  for  the  ith  contravariant 
component  of  equation  (1),  assuming  a  Newtonian  fluid,  is  given  by: 


p  u*7  y>7.  =  p  f  -  gU  P,.  +  u  v]  *  (\  +  V)  gU  v) 

J  T*  Jd 


ki 


2elJ*  v»j  +  gLZ  \ti 


(no  summation  on  i) 


(2) 


In  equation  (2),  e is  the  contravariant  form  of  the  deformation  tensor  and 
is  defined  by: 


h  (gzl  vL  *  g°'j  v*J 
*  c 


(3) 


The  quantity  y  is  the  fluid  viscosity  and  is  related  to  the  vatue  A  by 
Stokes'  hypothesis,  according  to  which: 


A  +  §  y  = 


0 


(4) 


Substitution  of  equation  (4)  into  equation  (2)  and  pre-multi plying  by  g ?. 


results  in: 


p  4*  **  =  p  Ai  f  -  jr  p*i +  v  gti  g"°  vhj + 


U  3  k  n  %  i  o  2  1k 

3  IT  V>ki  *  *Bii  e  v‘j  ~lTv»k  v*i 


Expressions  for  v,  .  and  u,  . .  are  needed  to  put  equation  (5)  into  a  more  useful 
0  00 

form.  The  quantity  v]  .  is  defined  as  the  jth  covariant  derivative  of  the  ith 

0 

contravariant  component  of  the  vector  velocity,  and  is  given  by  ^ 


12,  ♦  {  i .}  vm 
ZaP 


The  expression  {  1 .}  represents  the  second  kind  of  Christoffel  symbol  and  is 
m  j 

not  a  tensor  but  a  function  of  the  metric  tensor,  gtJ  ;  it  arises  due  to  the 

(2)  i 

variation  of  coordinate  base  vectors  •  Similarly,  the  form  v ,  ..  can  be 
,  .  3d 

(3) 

shown  to  be' 

4  ■  'ftps*  $***0  A  * 


•a  <A>  £  *  W  + 


Substitution  of  equations  (6)  and  (7)  into  equation  (5),  followed  by 

* 

considerable  algebraic  manipulation,  and  expressing  the  final  result  in  terms 


Details  of  this  derivation  are  available  in  reference  3 
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of  physical  components,  yields: 


-tf  7 h[9* p  v(i}  \  ~h  f1* gH  \  ^t\  =  sF(i)  *  sp(i)  + 

g*  te  1  J  J  g*  ZxJ  l  *2  hxO  )  F  P 

(8) 


h2.  dx° 

V 


0  J  g1 

SD(i)  -  SJi)  +  S  (i)  +  S  (i)  +  S  (i) 
u  c  P  V  PP 


where: 


Sp(i)  =  p  fd) 


o  _  1  3P 

p  "  tTTT 

V  3x 


(no  sum  on  i) 


and: 


(a)  Terms  arising  in  non-rectangular  geometries  only 
Term  Formula 


Sc(i)  =  p 


dh. 
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(10) 


(ID 


(12) 


I 


L 


i 


5 


(b)  Terms  arising  in  all  geometries  for  variable  density  and/or  viscosity 


Term 


Formula 


SJi) 

P 


-  V  J_ 

-  3h. 


v(k) 


3x 


3x 


■f  {  \\ 
k‘ 


k  >  v(m) 


m 


(13) 


svd) 


-  h . 


3y 


foh'.1  v(j) 


K 


♦  ~jl 
ln  i1  h 


Z  3^ 

fs/r.1  vd) 


3x 


3a^ 


■f  {  \} 


m  c‘  h 
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■m  J 


1  *v(i) 


h.  h2.  7>3? 


(14) 
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_  2  1  3y 
2  hi 


v(j)  .  ,  , 

^  m  **  hm 


(15) 


Equation  (8)  is  the  physical  component  form  of  the  momentum  equation  in 
general  curvilinear  orthogonal  coordinates  and  corresponds  to  equation  (4)  in 
reference  1.  It  has  been  formulated  here  such  that  effects  due  to  variation 
in  physical  properties  appear  in  explicitly  defined  terms,  SJi),  SJi) 

SpJi).  These  terms  represent  additional  contributions  to  the  momentum  balance 
and  are  volume  integrated  in  the  same  manner  as  SJi)  and  SJi)  in  reference  1 
to  obtain  curvilinear  coordinate  forms  of  the  difference  equations  in  variable 
property  flows. 

Specific  forms  for  the  source  terms  in  a  curvilinear  system  of  interest 
are  readily  derived  from  the  tabulated  expressions  simply  from  a  knowledge  of 
the  metric  tensor  for  that  curvilinear  system. 
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Notation 

Jj 


f 

fli) 

f 

H 


>3 


sc(i) 

sD(i) 

sFm 

spu) 

sva) 

SQ(i) 

SnJi) 

PM 


v 

V  (i) 

n 


deformation  tensor  in  contravariant  form  (i,j  =  1,  2,  3) 
component  of  body  force  in  tensor  notation  [i  =  1,  2,  3) 
physical  component  of  body  force  (i  =  1,  2,  3) 
body  force  field  in  vector  form 

L 

(=  (h  h  h  )*)  square  root  of  the  metric  tensor  determinant 

loo 

metric  tensor  in  contravariant  form  (i,j  =  1,  2,  3) 

scale  factor  in  curvilinear  orthogonal  coordinates  (£  =  1,  2,  3) 

pressure 

convection  source  term  in  momentum  equation  (t  8  1,  2,  3) 

diffusion  source  term  in  momentum  equation  (£  =  1,  2,  3) 

body  force  source  term  in  momentum  equation  (i  -  1,  2,  3) 

pressure  source  term  in  momentum  equation  (i  =  1,  2,  3) 

source  term  due  to  variable  viscosity  in  momentum  equation [i  =  1,  2,  3) 

source  term  due  to  variable  density  in  momentum  equation (f  =  1.  2,  3) 
source  term  due  to  variable  density  and  viscosity  in  momentum(i  =  1.  2,  3) 
stress  tensor  in  vector  notation 

velocity  component  (contravariant)  in  tensor  notation  C2  :=/.*»$) 
physical  component  of  velocity  (i  =  1,  2,  3) 
velocity  in  vector  notation 

coordinate  component  (contravariant)  in  tensor  notation  ( ^  ) 
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constant  related  to  p  by  Stokes'  Hypothesis  in  equation  (4) 

viscosity 

mass  density 

gradient  operator  (Hamilton's  operator) 
divergence  operator 

Christoff el  symbol  of  second  kind 
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